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ABSTRACT: Eighteen active substances, including 17 organosulfur 
compounds found in garlic essential oil (T), were identified by GC—MS 
analysis. For the first time, using the molecular docking technique, we report 
the inhibitory effect of the considered compounds on the host receptor 
angiotensin-converting enzyme 2 (ACE2) protein in the human body that 
leads to a crucial foundation about coronavirus resistance of individual 
compounds on the main protease (PDB6LU7) protein of SARS-CoV-2. The 
results show that the 17 organosulfur compounds, accounting for 99.4% 
contents of the garlic essential oil, have strong interactions with the amino 
acids of the ACE2 protein and the main protease PDB6LU7 of SARS-CoV- 

2. The strongest anticoronavirus activity is expressed in allyl disulfide and 

allyl trisulfide, which account for the highest content in the garlic essential 

oil (51.3%). Interestingly, docking results indicate the synergistic interactions of the 17 substances, which exhibit good inhibition of 

the ACE2 and PDB6LU7 proteins. The results suggest that the garlic essential oil is a valuable natural antivirus source, which 

contributes to preventing the invasion of coronavirus into the human body. 


1. INTRODUCTION 

Garlic (Allium sativum L.) (cf. Figure l) is considered as an 
important herb thanks to its variety of uses, including either a 
common spice for family meals or a popular component in 
folk-medicine prescriptions. 1 ' 2 For thousands of years, garlic 



Figure 1. Picture of garlic (A. sativum L.). 


has been used as a medication for common colds, influenza, 
and other kinds of infections. 3 ' 4 Recent pharmacological 
studies indicate that essential oil of garlic is an exceptional 
source of organosulfur compounds, possessing strong anti¬ 
oxidant, antibacterial, antifungal, anticancer, and antimicrobial 
properties. The oil is also proven to be conducive to 
hypoglycemia, hypotension, antithrombotic, immunomodula¬ 
tory, and prebiotic therapy. Besides, allicin is a typical reactive 
sulfur species found in the essential oil. 5 

Recently, many people have been infected with a novel 
coronavirus (SARS-CoV-2), and the death toll has reached 
thousands and been increasing day by day, which is a major 
problem in the world. 6 Therefore, the demand to seek for 
natural and safe medicines to prevent coronavirus is of great 
concern for all scientists around the world. With abundant 
medicinal resources in Vietnam and the specific healing 
properties of garlic/ we herein recommend a solution for the 
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prevention and treatment of the Coronavirus disease 2019 
(COVID-19) using the garlic essential oil. 

The fact is that coronaviruses belong to a large family of 
viruses that often cause the common cold in humans. Middle 
East respiratory syndrome (MERS), severe acute respiratory 
syndrome (SARS), and lately SARS-CoV-2 are the more severe 
symptoms caused by the coronavirus family. 6 ' 8 SARS-CoV-2 is 
a new strain that has been unprecedentedly found in humans. 
Angiotensin-converting enzyme 2 (ACE2) is an integral 
membrane glycoprotein that is known for the highest 
expression in most tissues such as kidneys, endothelium, 
lungs, and heart. 9, J The ACE2 protein is the same functional 
host-cell receptor shared by SARS-CoV-2 and SARS. 8,11,1 z The 
structural database of ACE2 can be referenced from 
UniProtKB. 1 ^ Therefore, besides inhibiting SARS-CoV-2, the 
inhibition of the ACE2 protein is absolutely necessary to 
reduce the operability of the host receptor of SARS-CoV-2. If 
the ACE2 protein is inhibited, it suggests that coronavirus is 
prevented and treated. 12 

As mentioned above, the high organosulfur compounds in 
garlic essential oil are expected to have strong interactions with 
the amino acids of the ACE2 protein. In this study, the idea is 
to determine the inhibition capacity of ligands in garlic 
essential oil not only to the ACE2 protein (host receptor for 
SARS-CoV-2) but also directly to the PDB6LU7 protein (main 
protease of SARS-CoV-2). 14 The structure of the PDB6LU7 
protein of SARS-CoV-2 has been determined recently by 
Worldwide Protein Data Bank. 14 It shows that the two proteins 
selected from Uniprot and Worldwide Protein Data Bank have 
three-dimensional structures, as shown in Scheme 1 with the 
DOI: 10.2210/pdbACE2/pdb presented for the ACE2 protein 
and the DOI: 10.2210/pdb6LU7/pdb shown for the 
PDB6LU7 protein in SARS-CoV-2. 11-14 


Scheme 1. (A) Angiotensin-Converting Enzyme 2 (ACE2) 
in the Human Body and (B) PDB6LU7 Protein in SARS- 
CoV-2 Main Protease 



In this study, garlic essential oil was extracted from 
commercial garlic and collected in Hue, Vietnam, by steam 
distillation. The composition of the essential oil was identified 
by GC—MS. The fact is that the study of HIV-1 resistance to 
reverse transcriptase inhibitors was reported by Tarasova et 
al. 15 The inhibition abilities of the garlic essential oil toward 
the ACE2 and PDB6LU7 proteins were determined using the 
molecular docking technique to investigate the interactions of 
ligands in the garlic essential oil with the ACE2 and PDB6LU7 
proteins. The docking results (i.e., docking score energy (DS), 
root-mean-square deviation (RMSD), distances, and types of 
interactions) indicate the inhibitory effects of the compounds 


in the garlic essential oil on the two proteins. To the best of 
our knowledge, this is the first report on the inhibitory effects 
of the considered garlic compounds on the ACE2 protein, 
which is a crucial foundation about SARS-CoV-2 resistance of 
those compounds on the main protease (PDB6LU7) protein 
of SARS-CoV-2 using docking simulations. The results of this 
study show that natural garlic essential oil is considered as a 
valuable resource recommended for preventing SARS-CoV-2 
invasion into the human body. 

2. RESULTS AND DISCUSSION 

2.1. Composition of Garlic Essential Oil. The density 
and refractive index of the garlic essential oil (A. sativum L.) 
are 1.019 g*mL _1 and 1.467, respectively. The results of the 
analysis of physicochemical properties in this study agree well 
with the results (1.025—1.029 g-mL _1 ; 1.464—1.469) reported 
by Boukeria et al. 16 

The qualitative and quantitative composition of the garlic 
essential oil was determined by GC—MS analysis, and the data 
are recorded in Figure 2 and Table 1, where the compounds 
are listed in order of their elution. A total of 18 compounds 
were identified in the garlic essential oil, covering more than 
96.6% of the GC—MS profiles. The main constituents in the 
garlic essential oil were allyl disulfide (28.4%), allyl trisulfide 
(22.8%), allyl (E)-l-propenyl disulfide (8.2%), allyl methyl 
trisulfide (6.7%), and diallyl tetrasulfide (6.5%). These results 
had differences in concentrations of chemical composition 
compared to those in the previous studies. 3,17-19 Li et al. 
indicated that major essential oil components of garlic from 
Suzhou City in China were 3-vinyl-4H-l,2-dithiin (31.9%), 
diallyl trisulfide (13.3%), diallyl sulfide (2.2%), diallyl disulfide 
(6.9%), propyl allyl disulfide (13.9%), and dimethyl disulfide 
(7.1%)/ On the contrary, diallyl trisulfide (37.3—45.9%), 
diallyl disulfide (17.5—35.6%), and methyl allyl trisulfide (7.7— 
10.4%) were the major components of the garlic essential oil 
from the state of Ariana (Tunisia). 19 The main essential oil 
components obtained from the east of Algeria were diallyl 
disulfide (33.4—38.4%) and diallyl trisulfide (23.35—30.0%). 20 
It has been known that the difference in the chemical 
composition of essential oils can be attributed to the 
importance of geo-ecological factors in the production of 
metabolites of plants. 3,17- 9 

Although the bioactivities of garlic essential oil have not 
been investigated in this study, they were documented in 
previous studies. 21,22 According to Rattanachaikunsopon 
Pongsak and Parichat Phumkhachorn, allyl disulfide (Tl) 
was also able to inhibit Escherichia coli 0157:H7 in a food 
model. 25 Moreover, Wodek et al. indicated that allyl disulfide 
can constitute a source of sulfane sulfur for liver; therefore, it 
can be used for cyanide detoxification in mouse tissues. 22 
Methyl allyl disulfide and diallyl trisulfide were inhibitory 
agents against Sitophilus zeamais Motschulsky and Triholium 
castaneum (Herbst) for contact toxicity, fumigant toxicity, and 
antifeedant activity. Besides, according to Seki et al., the 
anticancer effect of diallyl trisulfide on cancer cells HCT-15 
and DLD-1 was reported. 24 The biological activities of the 
essential oil depend on the number of sulfur atoms, and the 
higher the number of sulfur atoms, the stronger the biological 
activities. 1 ' In this study, the high amount of sulfur (95.63%) 
in the garlic essential oil suggests a further application of this 
essential oil in the medicine and pharmaceutical industry. 

2.2. Docking Simulation. To the best of our knowledge, 
there are no results related to SARS-CoV-2 resistance through 
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GC-MS for composition of garlic essential oil 



Figure 2. GC—MS chromatogram of garlic essential oil. 


theoretical and experimental studies of compounds in essential 
oils. 

2.2.1. Docking Simulation Results of Compounds in 
Garlic Essential Oil into ACE2 Protein. This is our first 
successful demonstration of the model describing the docking 
molecules of 17 out of 18 compounds in the garlic essential oil 
into the complex structures of the ACE2 protein. Docking 
simulations of the interactions between compounds Tl—T17 
and the ACE2 protein in the human body are presented as Tl- 
ACE2 to T17-ACE2 in the Supporting Information. We found 
that T18-ACE2 has the lowest content in the essential oil, and 
due to its bulky structure, the interaction with the ACE2 
protein is not easy; hence, we do not demonstrate its 
simulation in this study. The results of DS energy and 
RMSD between compounds in garlic essential oil and proteins 
with various interactions such as hydrogen bonds, cation—# 
bonds, 71—71 bonds, and ionic interactions, as well as the 
interaction distance between amino acids and the active sites of 
compounds are presented in Table 2, Figures 3, and SI 
(Supporting Information). There are also van der Waals 
interactions between compounds with other amino acids of the 
ACE2 protein, but since they are weak interactions, they are 
not identified in this study. 

The model describing the docking molecules indicated that 
all of the 17 compounds in the garlic essential oil from Tl- 
ACE2 to T17-ACE2 have the binding ability toward the areas 
affected by the ACE2 protein. The five presented representa¬ 
tives T5-ACE2, T11-ACE2, T1-ACE2, T2-ACE2, and T4- 
ACE2 have good DS energies of —14.06, —14.01, —12.84, 
— 12.76, and —12.5 kcahmol -1 , respectively (cf. Table 2 and 
Figure 3). It can be seen that T5-ACE2 interacts best with 
regions of the ACE2 protein, and the same trend is found for 
T11-ACE2. The priority order of the interactions of 
compounds in the garlic essential oil was as follows: T5 = 
Til > Tl = T2 > T4 > T8 > T9 > T12 > T13 > T14 > T15 > 
T3 > T7 > T10 > T16 > T17 > T6. 

We found that the sulfur compounds had strong interactions 
with amino acids in the ACE2 protein. It is revealed from the 
results of the model describing the docking molecules that the 
relation in terms of structure and inhibition ability of the 17 
substances in the garlic essential oil toward the ACE2 protein 
was as follows: the compounds had very strong interactions 
with the amino acids Pro 565, Trp 566, Ala 396, Gin 102, Gin 
101, Glu 208, Gly 205, Gin 98, Asn 210, Lys 94, Lys 562, Val 
209, and Ser 563. The interactions mainly took place with the 


molecule groups containing sulfur in their compounds. 
Interestingly, T10-ACE2 and T13-ACE2 mainly interacted 
with oxygen in the compounds. The total content of the 17 
compounds in the garlic essential oil resistant to the ACE2 
protein was 99.4%. This proves that the entire essential oil is 
appropriate for the prevention and treatment of pneumonia 
caused by SARS-CoV-2. Thus, finding active sites of 
compounds in the essential oil to inhibit the ACE2 protein 
is of great significance to the orientation of using the garlic 
essential oil in the prevention and treatment of SARS-CoV-2 in 
specific and other viruses causing flu or pneumonia in general. 

2.2.2. Docking Results of Compounds in Garlic Essential 
Oil into the PDB6LU7 Protein of SARS-CoV-2. The model 
describing the docking molecules of 18 substances in the garlic 
essential oil into the complex structure of the PDB6LU7 
protein was constructed using MOE 2015.10 software. The 
docking was successful in 17 compounds of Tl-SARS-CoV-2 
to T17-SARS-CoV-2. The DS energy and RMSD values 
between ligands and proteins shown through hydrogen bonds, 
cation—# bonds, #—# bonds, and ionic interactions, as well as 
the interaction distance between amino acids and the active 
sites in the title molecules are presented in Table 3 and Figures 
4 and S2 (Supporting Information). 

The anti-SARS-CoV-2 activity of the garlic essential oil 
composition into the PDB6LU7 protein of SARS-CoV-2 was 
found in the following order: T2 = Tl > T5 > T4 > Til > 
T15 > T8 > T16 > T9 > T12 > T13 > T3 > T6 > T7 > T10 > 
T14 > T17. The compounds Tl and T2 indicated the highest 
content in the garlic essential oil (51.3%) and gave the best 
anti-SARS-CoV-2 activity, followed by the three compounds 
T5, T4, and Til. These five substances account for 65.83% 
composition of the garlic essential oil (Figure 4). The results 
demonstrated that these 17 compounds have a stronger 
inhibition on PDB6LU7 of SARS-CoV-2 than that on the 
ACE2 protein in the human body. The DS energy was 
negative, in the range of —15.32 to —11.68 kcahmol -1 . 
Interactions with amino acid included Pro 565, Gin 102, Glu 
208, Asn 210, Gly 205, Gin 98, Trp 566, Lys 94, Val 209, Gin 
101, Asp 206, Asn 103, Ser 563, Ala 396, and Lys 562. The 
total content of 17 out of 18 compounds in the garlic essential 
oil exhibiting inhibition ability toward the PDB6LU7 protein 
of SARS-CoV-2 is 99.4%. The docking simulation results are 
excellent evidence of the anti-SARS-CoV-2 activity of the garlic 
essential oil. 
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Table 1. Identification of the Bioactive Compounds in Garlic Essential Oil 


No. 


10 


11 


Compound 


Formula Structure Symbol 


1 Diacetonalcohol 

2 Allyl sulfide 

3 Methyl allyl disulfide 

4 1,2-Dithiole 

5 Allyl disulfide 

6 Allyl (Z)-l-propenyl disulfide 

7 Allyl (E)-l-propenyl disulfide 

8 Allyl methyl trisulfide 


3-Vinyl- 1 ,2-dithiacyclohex-4-ene C 6 H 8 S 2 


2-Vinyl-4H-1,3-dithiine 


Carvone 


12 Allyl trisulfide 

13 Trisulfide, 2-propenyl propyl 


c 6 h 8 s 2 

Ci 0 H 14 O 

c 6 h 10 s 3 

c 6 h 12 s 2 


14 Trisulfide, (lZ)-l-propenyl 2-propenyl C 6 Hi 0 S 3 

15 Trisulfide, (IE)-1-propenyl 2-propenyl C 6 Hi 0 S 3 

16 1-propenyl methyl disulfide C 4 H 8 S 2 

17 Diallyl tetrasulfide C 6 Hi 0 S 3 


C 6 H 12 0 2 HO. 

c 6 h 10 s 
c 4 h 8 s 2 

C 3 H 4 S 2 
C 6 Hi 0 S 2 
C 6 HioS 2 
C 6 Hi 0 S 2 

C 4 H 8 S 3 /S'; 


's 




'S's'S' 


3 'S' S ' 


3 ' S ' S -^ 


18 


Cyclic octaatomic sulfur 




,S-S X 

s s 
s N ,s 
s-s 


T13 

T15 

T12 

T6 

T1 

T7 

T3 

T4 

T9 

T8 

T10 

T2 

Til 

T17 

T14 

T16 

T5 

T18 


Percentage (%) 

0.99 

0.97 

1.14 

5.21 

28.44 

3.19 

8.18 

6.74 

2.93 

3.11 

2.10 

22.81 

1.33 

0.64 

0.99 

0.68 

6.51 

0.60 


The simultaneous interaction of 17 out of 18 compounds in 
the garlic essential oil with the ACE2 protein and the 
PDB6LU7 protein is described in Figure 5. Hydrogen bonds, 
cation—;r bonds, n—n bonds, and ionic interactions of 17 
compounds with the amino acids of the ACE2 protein 
included Trp 566, Asn 103, Gly 205, Lys 94, Asp 206, Glu 
208, Val 209, Pro 565, Asn 210, Ser 563, Gin 98, Gin 102, Gin 
101, Ala 396, Lys 562, and Ala 396, and those with the amino 
acids of the PDB6LU7 protein included Gly 143, Glu 166, Asn 
142, Leu 141, Cys 145, Ser 144, His 163, and Met 165. Based 
on the percentage of each compound in the garlic essential oil 


(shown by GC—MS) and docking score (DS) energy of each 
compound, the average DS energies of these 17 compounds 
with the ACE2 protein and the PDB6LU7 protein were 
calculated as —11311 and —13 871 kcal-mol -1 , respectively. 
This shows an excellent result to confirm the effects of the 
garlic essential oil on the virus host receptor (ACE2) inhibition 
and SARS-CoV-2 resistance. Therefore, it is possible to use 
each compound in garlic essential oil or the whole essential oil 
system to act simultaneously on the ACE2 protein and the 
PDB6LU7 protein. This result opens the direction of research 
and application of the essential oils in general, garlic essential 
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Table 2. Docking Simulation Results with Docking Score Energy (DS) and Root-Mean-Square Deviation (RMSD) between 
Compounds in Garlic Essential Oil and the ACE2 Protein 


T5 = Til > T1 = T2 > T4 > T8 > T9 > T12 > T13 > T14 > T15 > T3 > T7 > T10 > T16 > T17 > T6 


compound 

symbol (compound- 
protein) 

DS ( 

kcal-mol -1 ) 

RMSD 

(A) 

interaction with amino acid 

allyl disulfide 

T1-ACE2 

-12.84 

1.46 

Pro 565 (2.85 A), Trp 566 (2.87 A), Ala 396 (3.14 A) 

allyl trisulfide 

T2-ACE2 

-12.76 

1.47 

Gin 102 (2.89 A; 3.11 A), Asn 103 (3.45 A) 

allyl (E)-l-propenyl disulfide 

T3-ACE2 

-9.07 

0.66 

Glu 208 (3.06 A), Gly 205 (3.70 A) 

allyl methyl trisulfide 

T4-ACE2 

-12.50 

1.56 

Asn 210 (3.13 A; 3.05 A), Lys 94 (3.22 A) 

diallyl tetrasulfide 

T5-ACE2 

-14.06 

1.23 

Gly 205 (3.43 A), Asp 206 (2.95 A), Lys 562 (2.92 A), Ala 396 (3.82 A) 

1,2-dithiole 

T6-ACE2 

-7.89 

3.09 

Asn 210 (3.29 A) 

allyl (Z)-l-propenyl disulfide 

T7-ACE2 

-9.04 

1.58 

Gin 98 (4.14 A), Glu 208 (3.24 A) 

2-vinyl-4H- 1,3-dithiine 

T8-ACE2 

-11.83 

0.62 

Gin 98 (3.10 A), Val 209 (3.56 A), Asn 210 (3.93 A) 

3-vinyl- l,2-dithiacyclohex-4-ene 

T9-ACE2 

-10.57 

1.19 

Trp 566 (2.78 A), Pro 565 (3.95 A), Glu 208 (3.12 A) 

carvone 

T10-ACE2 

-8.58 

1.53 

Lys 94 (2.13 A), Gin 98 (1.72 A) 

trisulfide, 2-prop enyl propyl 

T11-ACE2 

-14.01 

1.85 

Tip 566 (3.51 A), Glu 208 (3.66 A), (3.05 A), Val 209 (3.25 A), Gin 98 
(3.52 A) 

methyl allyl disulfide 

T12-ACE2 

-10.32 

1.16 

Val 209 (2.31 A), Gin 98 (3.45 A), Lys 94 (3.15 A) 

diacetonalcohol 

T13-ACE2 

-9.71 

0.85 

Gin 101 (2.13 A), Asn 210 (1.97 A; 2.05 A) 

trisulfide, (lE)-l-propenyl 2- 
propenyl 

T14-ACE2 

-9.57 

0.68 

Asn 210 (3.20 A; 1.72 A), Ser 563 (3.73 A) 

allyl sulfide 

T15-ACE2 

-9.38 

1.48 

Asn 210 (2.31 A), Lys 94 (3.15 A), Gin 98 (3.45 A) 

1-propenyl methyl disulfide 

T16-ACE2 

-8.06 

0.88 

Asp 206 (2.81 A), Trp 566 (3.46 A) 

trisulfide, (1Z)-1-propenyl 2- 
propenyl 

T17-ACE2 

-8.06 

2.35 

Glu 208 (3.29 A), Gin 98 (3.64 A) 


Docking Study for Angiotensin-Converting Enzyme 2 (ACE2) 




T5-ACE2 


O polar — * sidechain acceptor O solvent residue arene-arene 

O acidic sidechain donor O metal complex @H arene-H 
O basic — * backbone acceptor solvent contact arene-cation 

O greasy •* backbone donor — metal/ion contact 
proximity bgand p. receptor 

contour exposure ^ exposure 



T11-ACE2 


T1-ACE2 


T2-ACE2 T4-ACE2 


Figure 3. (A) Native human angiotensin-converting enzyme 2 (ACE2) crystal structure. Docking simulation of the interaction between 
compounds T5, Til, Tl, T2, T4, and the ACE2 protein in the human body: (B) T5-ACE2, (C) T11-ACE2, (D) T1-ACE2, (E) T2-ACE2, and 
(F) T4-ACE2. The inhibitory effects of the compounds on the ACE2 protein are of the order: T5 « Til >T1 «T2> T4. 


oil in particular, in the prevention and treatment of SARS- 
CoV-2. 

3. CONCLUSIONS 

This study proposes a potential approach to the use of natural 
essential oils, in general, and garlic essential oil, in particular, to 
tackle the current pandemic SARS-CoV-2. The compounds in 
the garlic essential oil inhibit the ACE2 protein, leading the 
virus to lose the host receptor and attacking the PDB6LU7 


protein—the main protease of SARS-CoV-2—at the same 
time. This prevents protein maturation of the virus and the 
spread of infection. Docking simulation suggests the active 
binding site of most active compounds in garlic essential oil 
with the ACE2 protein and the PDB6LU7 protein. From the 
analysis of the docking data, it is revealed that 17 (Tl—T17) 
out of 18 compounds of the garlic essential oil are capable of 
inhibiting ACE2 and resisting SARS-CoV-2 and that the total 
content of these 17 compounds accounts for 99.4% 
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Table 3. Docking Simulation Results with Docking Score Energy (DS) and Root-Mean-Square Deviation (RMSD) between the 
Title Compounds and PDB6LU7 Protein of SARS-CoV-2 


T2 = Tl > T5 > T4 > Til > T15 > T8 > T16 > T9 > T12 > T13 > T3 > T6 > T7 > T10 > T14 > T17 

symbol (com] 
compound protein 

pound- DS ( RMSD 

j kcal-mol -1 ) (A) interaction with amino acid 


allyl disulfide 

Tl-SARS-CoV-2 

-15.32 

1.35 

allyl trisulfide 

T2-SARS-CoV-2 

-15.02 

0.66 

allyl (E)-1-propenyl disulfide 

T3-SARS-CoV-2 

-13.25 

1.49 

allyl methyl trisulfide 

T4-SARS-CoV-2 

-14.36 

1.37 

diallyl tetrasulfide 

T5-SARS-CoV-2 

-14.47 

0.94 

1,2-dithiole 

T6-SARS-CoV-2 

-13.21 

2.30 

allyl (Z)-l-propenyl disulfide 

T7-SARS-CoV-2 

-12.60 

1.72 

2-vinyl-4H- 1,3-dithiine 

T 8 -SARS-C 0 V -2 

-14.04 

4.35 

3 -vinyl-1,2- dithiacyclohex-4- 

T9-SARS-CoV-2 

-13.83 

1.05 

ene 




carvone 

T10-SARS-CoV-2 

-12.36 

1.62 

trisulfide, 2-propenyl propyl 

Tll-SARS-CoV-2 

-14.36 

1.37 

methyl allyl disulfide 

T12-SARS-CoV-2 

-13.56 

1.51 

diacetonalcohol 

T13-SARS-CoV-2 

-13.26 

1.41 

trisulfide, (IE)-1-propenyl 2- 
propenyl 

T14-SARS-CoV-2 

-12.00 

2.12 

allyl sulfide 

T15-SARS-CoV-2 

-14.24 

1.20 

1-propenyl methyl disulfide 

TI 6 -SARS-C 0 V -2 

-13.84 

0.91 

trisulfide, (1Z)-1-propenyl 
2-propenyl 

T17-SARS-CoV-2 

-11.68 

3.60 


Gly 143 (2.59 A), Asn 142 (2.92 A), Leu 141 (2.94 A), Cys 145 (2.98 A), Ser 144 
(3.15 A), His 163 (2.95 A) 

Gly 143 (3.93 A; 2.79 A), Asn 142 (2.87 A; 2.67 A), Cys 145 (3.25 A; 3.63 A), Ser 
144 (3.06 A) 

Leu 141 (2.98 A), Ser 144 (3.07 A) 

Gly 143 (3.69 A), Asn 142 (2.81 A), Cys 145 (3.61 A; 3.12 A), Ser 144 (2.96 A) 

Gly 143 (2.88 A; 2.98 A), Asn 142 (2.10 A), Cys 145 (3.94 A; 3.56 A), Ser 144 
(3.14 A) 

Asn 142 (3.75 A), Cys 145 (3.75 A) 

Gly 143 (3.32 A), Leu 141 (3.62 A) 

Gly 143 (2.73 A), Asn 142 (3.82 A), Cys 145 (3.11 A), Ser 144 (3.02 A) 

Glu 166 (2.86 A), Met 165 (2.88 A), Cis 145 (2.13 A) 

His 163 (2.15 A) 

Gly 143 (2.62 A; 2.81 A), Asn 142 (2.85 A), Cys 145 (3.15 A; 3.50 A), Ser 144 
(2.87 A), His 163 (2.97 A) 

Cys 145 (2.43 A), Ser 144 (2.82 A), His 163 (2.65 A) 

Cys 145 (2.93 A), Ser 144 (2.26 A) 

Leu 141 (3.32 A), Ser 144 (2.39 A) 

Gly 143 (2.75 A), Asn 142 (2.75 A), Leu 141 (3.56 A), Cys 145 (3.47 A), Ser 144 
(3.32 A) 

Gly 143 (2.89 A), Asn 142 (3.88 A), Cys 145 (2.83 A), Ser 144 (2.80 A) 

Leu 141 (3.38 A) 
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Figure 4. (A) Crystal structure of the virus main protease in the complex (PDB6LU7). Docking simulation of the interaction between compounds 
Tl, T2, T4, T5, and Til, and the PDB6LU7 protein of SARS-CoV-2: (B) T1 -SARS-CoV-2, (C) T2-SARS-CoV-2, (D) T4-SARS-CoV-2, (E) 
T5-SARS-CoV-2, and (F) Tll-SARS-CoV-2. The inhibitory effects of the compounds on the PDB6LU7 protein of SARS-CoV-2 were of the 
order: Tl « T2 > T5 > T4 > Til. 


composition of the garlic essential oil. The docking score (DS) 
energy of compounds for the ACE2 protein ranges from 
— 14.06 to —7.89 kcal-mol -1 , and the DS energy of compounds 
for the PDB6LU7 protein of SARS-CoV-2 ranges from —15.32 


to —11.68 kcal-mol \ The order of active compounds 
inhibiting the ACE2 protein is as follows: T5 = T11>T1 = 
T2 > T4 > T8 > T9 > T12 > T13 > T14 > T15 > T3 > T7 > 
T10 > T16 > T17 > T6. Meanwhile, the order of active 
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ACE2 SARS-CoV-2 



Figure 5. 3D-docking results with simultaneous interaction of all 17 substances (17T = Tl—T17) in the garlic essential oil with (A) ACE2 protein 
and (B) PDB6LU7 protein of SARS-CoV-2. 


compounds resisting SARS-CoV-2 is T2 = Tl > T5 > T4 > 
Til > T15 > T8 > T16 > T9 > T12 > T13 > T3 > T6 > T7 > 
T10 > T14 > T17. The synergistic interactions of 17 
substances of the garlic essential oil exhibited good inhibition 
on the ACE2 protein (host receptor of the virus) and the 
PDB6LU7 protein of the virus. This study opens the door 
toward the use of the garlic essential oil in discovering and 
treating SARS-CoV-2 to prevent the current pandemic. 

4. MATERIALS AND METHODS 

4.1. Experiment. 4.1.1. Sample Collection. Garlic (A. 
sativum L.) was purchased from a local market (Thua Thien 
Hue province in Central Vietnam). Then, it was botanically 
identified, and a voucher specimen was deposited at the 
Department of Biology, University of Sciences, Hue University. 
During the experiment, all of the samples were maintained in 
our lab at appropriate conditions (e.g., 25 °C in the dark). 

4.1.2. Chemicals and Equipment. All analytical-grade 
chemicals obtained from Sigma-Aldrich Co. were used without 
any further purification. The major equipment includes a 
hydrodistillation apparatus (Merck Specialities Pvt. Ltd., 
India), a polarimeter (Reichert Cat #14003000), a Jasco V- 
630 spectrophotometer (Japan), and a gas chromatography- 
mass spectrometry (GC—MS) device (Aligent GC 7890B-MS 
5975C). 

4.1.3. Isolation of Essential Oil. A hydrodistillation 
Clevenger apparatus was used to extract the garlic essential 
oil. Peeled garlic (250 g) was minced and then put into a 1 L 
flask, in which 400 mL of distilled water was added. The 
essential oil was obtained by water distillation at 100 °C for 90 
min, according to the Vietnamese Pharmacopoeia. 25 A 
sterilized glass vial was used to collect the essential oil. The 
essential oil was dried by anhydrous Na 2 S0 4 before being 
stored at 4 °C for further analysis. The extraction process was 
repeated three times to check the repeatability of the analytical 
procedure. 25 

4.1.4. Refractive Index and Density of the Essential Oil. 
Refractive index and density are two common physical 
parameters for checking essential oil composition and purity. 
The refractive index of essential oils was determined using a 
polarimeter, and the density of the essential oil was determined 
using a pycnometer flask according to the Vietnamese 
Pharmacopoeia 25 and Boukeria et al. 16 

4.1.5. Gas Chromatography-Mass Spectrometry (GC- 
MS) Analysis. An Agilent GC 7890B-MS 5975C instrument 
combined with an HP-5MS column (30 m X 250 pm X 0.25 
pm) was used to identify the chemical composition of the 


garlic essential oil. The utilized carrier gas was helium at 
constant pressure (13 psi). The essential oil (l pL) was 
injected into the GC with a split ratio of 20:1. The injection 
temperature was 250 °C. The column temperature program 
was applied as follows: start at 70 °C and then increased at the 
rate of 10 °C every minute until it reached 280 °C. After the 
analytes were separated on a capillary column, they passed 
through the ionization zone in the MS source (ionization 
energy: 70 eV; interface temperature: 280 °C; MS temper¬ 
ature: 230 °C; quadrupole temperature: 150 °C), and the 
neutral molecules were ionized yielding specific mass/charge 
(ra/z) ratios. C8—C30 Alkanes Calibration Standard (Sigma- 
Aldrich) was applied to identify unknown compounds through 
their retention indices by comparing their mass spectra to 
those contained in the NIST02 database. The concentration of 
each analyte will be calculated based on its peak area in the 
chromatogram. 

4.2. Molecular Docking Simulation. The molecular 
docking technique is used to predict and describe the 
interaction (including the binding energy and structural 
parameters) of compounds resistant to the ACE2 protein in 
the human body and the PDB6LU7 protein in SARS-CoV-2. 
The docking results reveal that there are potential compounds 
in the essential oil to prevent and treat SARS-CoV-2. The 
molecular docking modeling consists of five essential 
steps as follows: 

Step 1: Selection of proteins 

• Selection of proteins in PDB: The biological targets are 
the ACE2 protein and the PDB6LU7 protein presented 
at Uniprot 1 ^ and Worldwide Protein Data Bank. 14 

• Determination of bonding position: The protein action 
areas were determined based on the ligand positions 
within a radius of 4.5 A and the presence of important 
amino acids. Water molecules were removed, and the 
structures of amino acids were checked before 
reestablishing the enzyme action areas. 

Step 2: Preparation of proteins 

• Creating 3D structures: The 2D chemical structure (flat 
structure) of the compounds was automatically con¬ 
verted to the 3D chemical structure (three-dimensional 
structure) by ChemBioOffice 2018 software. 

• Lowest-energy state: The 3D molecular structures of the 
compounds were then optimized using SYBYL-X 1.1 
software with the purpose to correct the mismatch 
values of bond lengths, bond angles, bending angles, and 
unusual nonbonding interactions due to the atoms in 
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different parts of the compounds occupying the same 
space. The energy minimization program used the Conj 
Grad method (conjugate gradient) to choose the stop 
point where the energy change is less than the cutoff 
value of 0.001 kcal-mol -1 with partial charges of 
Gasteiger—Huckel and the maximum number of 
iterations needed was 10 000. This method accumulates 
the information about the energy function from each 
iteration, which is suitable for small and large molecules. 
Molecular dynamics conduction was used to obtain the 
structures with the lowest energy via the principle of 
simulated annealing in Sybyl-X 1.1 software. The 
molecular structure was heated at a high temperature 
(700 K) for a certain time (1000 ps) so that the 
molecule rearranges its current state and then cooled to 
200 K for another 1000 ps to a steady state for 
computing the final configuration. The program 
automatically iterated five times to find the various 
configurations needed, and then the minimized energy 
was observed again to determine steric energy. In this 
way, the structural parameters of compounds with 
different energies were found, which were more durable 
than the original structures. 

Step 3: Redocking 

Redocking of protein-compound cocrystal structures: the 
redocking of protein—ligand complex cocrystal structures aims 
to assess the suitability of docking parameters. The process was 
carried out with three structures of compounds as follows: 

(1) Separation of compounds from homogenized complexes 
in proteins. 

(2) Separation of compounds from homogeneous complexes 
and repreparation using Sybyl-X 1.1 software. 

(3) Preparation of new compounds (structure drawing, 
structural optimization parameters on minimal energy, 
molecular dynamics calculations). 

Root-mean-square deviation (RMSD) values reflect the 
deviation of compounds' structures after docking compared to 
available structures in the crystal structure and comparing the 
interactions of compounds in the crystal structure after 
docking. The docking results are considered reliable when 
the RMSD value is <1.5 A, and the interactions between 
compounds and the initial enzyme are not significantly 
different. 

Step 4: Molecules docking into protein: docking inves¬ 
tigations 

• Docking of compounds into prepared proteins and 
compounds via conducting the docking process with 
MOE 2015.10 software with the following options: the 
method of placing compound fragments into the triangle 
matching; the maximum number of results for each 
iteration is 1000; the maximum number of results for 
each compound fragmentation is 200; retain the five best 
configurations of each compound in the bonding 
complex for further analysis. The best configuration is 
the one with the lowest docking score (DS) energy 
(kcal-mol -1 ). This score is the total energy consumed for 
the formation of bonding interactions between the 
compounds and the two proteins selected. 

Step 5: Docking results analysis 

• Evaluation of docking score (DS): Analysis of the 
interactions between the compounds and targeted 


proteins, and performance of interaction on 2D and 
3D planes using MOE 2015.10 software. Various 
interactions, such as van der Waals interactions, 
hydrogen bonds, cation—# bonds, n—n bonds, and 
ionic interactions, and the interaction distance between 
amino acids and the active sites of compounds are 
plotted. Van der Waals interactions are detected by 
contact with hydrophilic and hydrophobic surfaces 
between the compounds and the bonding point. 

■ ASSOCIATED CONTENT 

O Supporting Information 

The Supporting Information is available free of charge at 

https: / / pubs.acs.org/ doi /10.1021/ acsomega.0c00772. 

Docking simulation of the interaction between com¬ 
pounds and the ACE2 protein in human body; 
inhibitory effects of compounds on the ACE2 protein: 
T8 > T9 > T12 > T13 > T14 > T15 > T3 > T7 > T10 

> T16 > T17 > T6; docking simulation of the 
interaction between compounds and the PDB6LU7 
protein of SARS-CoV-2; and inhibitory effects of 
compounds on the PDB6LU7 protein: T15 > T8 > 
T16 > T9 > T12 > T13 > T3 > T6 > T7 > T10 > T14 

> T17 (PDF) 

■ AUTHOR INFORMATION 

Corresponding Authors 

Duong Tuan Quang — Department of Chemistry, University of 
Education, Hue University, Hue City 530000, Vietnam; 

© orcid.org/0000-0002-4926-0271; Email: dtquang@> 
hueuni.edu.vn 

Nguyen Thi Ai Nhung — Department of Chemistry, University 
of Sciences, Hue University, Hue City 530000, Vietnam; 

© orcid.org/0000-0002-5828-7898; Email: ntanhung(S> 
hueuni.edu.vn 

Authors 

Bui Thi Phuong Thuy — Faculty of Fundamental Science, Van 
Fang University, Ho Chi Minh City 700000, Viet Nam 
Tran Thi Ai My — Department of Chemistry, University of 
Sciences, Hue University, Hue City 530000, Vietnam 
Nguyen Thi Thanh Hai — Department of Chemistry, University 
of Sciences, Hue University, Hue City 530000, Vietnam 
Le Trung Hieu — Department of Chemistry, University of 
Sciences, Hue University, Hue City 530000, Vietnam 
Tran Thai Hoa — Department of Chemistry, University of 
Sciences, Hue University, Hue City 530000, Vietnam 
Huynh Thi Phuong Loan — Department of Chemistry, 
University of Sciences, Hue University, Hue City 530000, 
Vietnam; © orcid.org/0000-0002-3659-4571 
Nguyen Thanh Triet — Faculty of Traditional Medicine, 
University of Medicine and Pharmacy at Ho Chi Minh City, Ho 
Chi Minh City 700000, Vietnam 
Tran Thi Van Anh — Faculty of Pharmacy, University of 
Medicine and Pharmacy at Ho Chi Minh City, Ho Chi Minh 
City 700000, Vietnam 

Phan Tu Quy — Department of Natural Sciences & Technology, 
Toy Nguyen University, Buon Ma Thuot City 630000, Vietnam 
Pham Van Tat — Department of Environmental Engineering 
Hoa Sen University, Ho Chi Minh City 700000, Vietnam 


8319 


https://dx.doi.Org/10.1021 /acsomega.0c00772 
ACS Omega 2020, 5, 8312-8320 




ACS Omega 


http://pubs.acs.org/journal/acsodf 


Article 


Nguyen Van Hue — Faculty of Engineering and Food 
Technology, University of Agriculture and Forestry, Hue 
University, Hue City 530000\ Vietnam 
Nguyen Tien Trung — Laboratory of Computational Chemistry 
and Modeling Department of Chemistry, Quy Nhon University, 
Quy Nhon City 590000, Vietnam ; © orcid.org/0000-0001- 
5710-5776 

Vo Thanh Tung — University of Sciences, Hue University, Hue 
City 530000, Vietnam 

Lam K. Huynh — Department of Chemical Engineering 
International University, Ho Chi Minh City 700000, Vietnam-, 
Vietnam National University, Ho Chi Minh City 700000, 
Vietnam; © orcid.org/0000-0002-3683-2787 

Complete contact information is available at: 
https: / / pubs.acs.org/10.1021/ acsomega.0c00772 

Notes 

The authors declare no competing financial interest. 

■ ACKNOWLEDGMENTS 

This research was supported by the Vietnam National 
Foundation for Science and Technology Development 
(NAFOSTED). The authors would like to thank Prof. Jong 
Seung Kim (Korea University) for his valuable comments and 
discussion. Also, they highly appreciate Dr Mingle Li (Korea 
University) and Khue Lai (IU) for helping with figure 
preparation. 

■ REFERENCES 

(1) Satyal, P.; Craft, J. D.; Dosoky, N. S.; Setzer, W. N. The chemical 
compositions of the volatile oils of garlic ( Allium sativum) and wild 
garlic ( Allium vineale). Foods 2017, 6, 63. 

(2) Chung, L. Y. The antioxidant properties of garlic compounds: 
allyl cysteine, alliin, allicin, and allyl disulfide. /. Med. Food 2006, 9, 
205-213. 

(3) Najjaa, H.; Neffati, M.; Zouari, S.; Ammar, E. Essential oil 
composition and antibacterial activity of different extracts of Allium 
roseum L., a North African endemic species. C. R. Chim. 2007, 10, 
820-826. 

(4) Romeilah, R. M.; Fayed, S. A.; Mahmoud, G. I. Chemical 
compositions, antiviral and antioxidant activities of seven essential 
oils. J. Appl. Sci. Res. 2010, 6, 50—62. 

(5) Chekki, R. Z.; Snoussi, A.; Hamrouni, I.; Bouzouita, N. Chemical 
composition, antibacterial and antioxidant activities of Tunisian garlic 
( Allium sativum) essential oil and ethanol extract. Mediterr. J. Chem. 
2014, 3, 947-956. 

(6) Chen, J. Pathogenicity and transmissibility of 2019-nCoV—a 
quick overview and comparison with other emerging viruses. Microbes 
Infect. 2020, 22, 69-71. 

(7) Nguyen, T. C.; Chu, T. A. P.; Van Nguyen, H. In vitro 
evaluation of antibacterial activity of garlic Allium sativum against 
poultry pathogens and effect of garlic supplementation on ducking 
growth performance. Vietnam J. Sci. Technol. 2019, 57, 9. 

(8) Tipnis, S. R; Hooper, N. M.; Hyde, R; Karran, E.; Christie, G.; 
Turner, A. J. A human homolog of angiotensin-converting enzyme 
cloning and functional expression as a captopril-insensitive carbox- 
ypeptidase. J. Biol. Chem. 2000, 275, 33238—33243. 

(9) Tikellis, C.; Thomas, M. Angiotensin-converting enzyme 2 
(ACE2) is a key modulator of the renin angiotensin system in health 
and disease. Int. J. Pept. 2012, 2012, No. 256294. 

(10) Roberts, A.; Deming, D.; Paddock, C. D.; Cheng, A.; Yount, B.; 
Vogel, L.; Herman, B. D.; Sheahan, T.; Heise, M.; Genrich, G. L.; 
et al. A mouse-adapted SARS-coronavirus causes disease and 
mortality in BALB/c mice. PLoS Pathog. 2007, 3, No. e5. 

(11) Gralinski, L. E.; Menachery, V. D. Return of the Coronavirus: 
2019-nCoV. Viruses 2020, 12, 135. 


(12) Paraskevis, D.; Kostaki, E. G.; Magiorkinis, G.; 
Panayiotakopoulos, G.; Sourvinos, G.; Tsiodras, S. Full-genome 
evolutionary analysis of the novel corona virus (2019-nCoV) rejects 
the hypothesis of emergence as a result of a recent recombination 
event. Infect., Genet. Evol. 2020, 79, No. 104212. 

(13) Homo Sapiens (Human), https://www.uniprot.org/uniprot/ 
Q9BYF1. 

(14) The crystal structure of COVED-19 main protease in complex 
with an inhibitor N3. https://www.wwpdb.org/pdb?id=pdb_ 
000061u7. Initial deposition on 26 January 2020. 

(15) Tarasova, O.; Poroikov, V.; Veselovsky, A. Molecular docking 
studies of HIV-1 resistance to reverse transcriptase inhibitors: Mini¬ 
review. Molecules 2018, 23, 1233. 

(16) Boukeria, S.; Kadi, K.; Kalleb, R; Benbott, A.; Bendjedou, D.; 
Yahia, A. Phytochemical and physicochemical characterization of 
Allium sativum L. and Allium cepa L. Essential oils. J. Mater. Environ. 
Sci. 2016, 7, 2362-2368. 

(17) Li, R. ; Chen, W. C.; Wang, W. P.; Tian, W. Y.; Zhang, X. G. 
Extraction of essential oils from garlic ( Allium sativum) using ligarine 
as solvent and its immunity activity in gastric cancer rat. Med. Chem. 
Res. 2010, 19, 1092-1105. 

(18) El-Sayed, H. S.; Chizzola, R.; Ramadan, A. A.; Edris, A. E. 
Chemical composition and antimicrobial activity of garlic essential 
oils evaluated in organic solvent, emulsifying, and self-micro- 
emulsifying water based delivery systems. Food Chem. 2017, 221, 
196-204. 

(19) Dziri, S.; Casabianca, H.; Hanchi, B.; Hosni, K. Composition of 
garlic essential oil ( Allium sativum L.) as influenced by drying method. 
J. Essent. Oil Res. 2014, 26, 91—96. 

(20) Boubechiche, Z.; Chihib, N.-E.; Jama, C.; Hellal, A. 
Comparison of Volatile Compounds Profile and Antioxydant Activity 
of Allium sativum Essential Oils Extracted using Hydrodistillation, 
Ultrasound-Assisted and Sono-Hydrodistillation Processes. Indian J. 
Pharm. Educ. Res. 2017, 51, S281—S285. 

(21) Rattanachaikunsopon, P.j Phumkhachorn, P. Diallyl sulfide 
content and antimicrobial activity against food-borne pathogenic 
bacteria of chives ( Allium schoenoprasum). Biosci., Biotechnol., Biochem. 
2008, 72, 2987-2991. 

(22) W 3 odek, L. Alyl disulfide as donor and cyanide as acceptor of 
sulfane sulfur in the mouse tissues. Pharmacol. Rep. 2005, 57, 212— 
218. 

(23) Huang, Y.; Chen, S. X.; Ho, S. H. Bioactivities of methyl allyl 
disulfide and diallyl trisulfide from essential oil of garlic to two species 
of stored-product pests, Sitophilus zeamais (Coleoptera: Curculioni- 
dae) and Tribolium castaneum (Coleoptera: Tenebrionidae). /. Econ. 
Entomol. 2000, 93, 537—543. 

(24) Seki, T.; Hosono, T.; Hosono-Fukao, T.; Inada, K.; Tanaka, R; 
Ogihara, J.; Aiga, T. Anticancer effects of diallyl trisulfide derived 
from garlic. Asia Pac. J. Clin. Nutr. 2008, 17, 249—252. 

(25) Vietnamese Pharmacopoeia-, Medical Publishing House: Hanoi, 
Vietnam, 1997. 

(26) Babu, T. M. C.; Rajesh, S. S.; Bhaskar, B. V.; Devi, S.; 
Rammohan, A.; Sivaraman, T.; Rajendra, W. Molecular docking, 
molecular dynamics simulation, biological evaluation and 2D QSAR 
analysis of flavonoids from Syzygium alternifolium as potent anti- 
Helicobacter pylori agents. RSC Adv. 2017, 7, 18277—18292. 

(27) Ngo, T. D.; Tran, T. D.; Le, M. T.; Thai, K. M. Computational 
predictive models for P-glycoprotein inhibition of in-house chalcone 
derivatives and drug-bank compounds. Mol. Diversity 2016, 20, 945— 
961. 

(28) Thai, K. M.; Le, D. P.; Tran, T. D.; Le, M. T.; et al. 
Computational assay of Zanamivir binding affinity with original and 
mutant influenza neuraminidase 9 using molecular docking. J. Theor. 
Biol. 2015, 385, 31-39. 


8320 


https://dx.d 0 i. 0 rg/l 0.1021 /acsomega.0c00772 
ACS Omega 2020, 5, 8312-8320 




